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1 Introduction 



The lifetime of a free neutron is small by astrophysical standards and amounts to only 
about 900 s. Therefore, neutrons can affect significantly the plasma dynamics and kinetics 
only in such objects in which the decay of free neutrons is compensated by processes of 
neutron accumulation. In fact, there are two such processes: fission of helium (or nuclei 
of heavier elements) and proton-to-neutron conversion due to the weak or, at very high 
energies, electromagnetic interaction. Of course, we do not touch on the specific case of 
superdense matter inside neutron stars where neutrons are practically stable and their 
role is dominant in many respects. 

Very little attention was paid up to now to the role of free neutrons in the dynamics and 
radiation of nondegenerate astrophysical plasmas. This can be explained in part by the 
widely accepted view that unstable particles cannot influence significantly the majority of 
processes in a plasma except fast nuclear reactions. In this respect, we should especially 
emphasize the synthesis of primordial helium in the early Universe, in which free neutrons 
play the key role. However, this has little in common with problems of plasma physics 
itself. The same is true as far as the nucleosynthesis inside ordinary stars is concerned. 

The few exceptions from this tradition are the following. A number of authors (e.g., ^) 
considered neutrons created due to inelastic interaction of photons and protons (p + 7 — > 
n + TT"*") as a species, independent of the magnetic-field influence, which is capable of 
transferring energy in active galactic nuclei. In this case, the typical energies of particles 
exceed 10^^ eV, i.e., problems of this kind are related to the physics of cosmic rays and 
thus suffer ambiguities in formulation and models adopted. 

In [2, 3], advection-dominated accretion flows (ADAFs) in which ion temperature is 
very high were proposed as another class of sources where free neutrons could appear. 
In this case, neutrons are created due to fission of helium nuclei in a hot nonisothermal 
plasma [4, 5]. In the cited papers, the role of neutrons was only discussed in relation to 
generation of gamma-ray lines after neutron capture by heavy elements or generation of 
bremsstrahlung radiation due to proton-neutron collisions [Q, though the effect of neutron 
component on the plasma viscosity was considered, as well [6]. 

In this paper, we consider a wider class of phenomena related to the existence of free 
neutrons in plasmas created by compact objects. It will be shown that the list of objects 
in which neutrons appear in a natural way should be extended to include also relativistic 
jets in microquasars [7] and gamma-ray bursts. On the other hand, the conditions of the 
formation of the neutron component appear to be much less stringent than was believed. 

Let us recall the main properties of gamma-ray bursts (see also the reviews [8, 9]). The 
burst sources are distributed isotropically in the celestial sphere, have angular sizes less 
than the angular resolution of modern detectors, and produce (probably nonrecurrent) 
bursts of durations 0.1 — 100 s. The maximum of the spectrum falls in the gamma-ray 
band (about 200 keV), which explains the name of this phenomenon. Recent observa- 
tions of optical components of gamma-ray bursts [10-12] revealed that GRB sources are 
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located at distances comparable to the size of the Universe. The total energy release of a 
typical source amounts to 10^^—10^^ ergs. Consequently, the observed radiation should be 
generated in an ultrarelativistic flow of matter moving with a Lorentz factor F > 100 |T3[| , 
which allows one to avoid strong two-photon absorption. 

When discussing the physics of a two-component relativistic flow (flreball), we will 
focus on issues not considered in detail in our previous papers [14, 15]. In particular, we 
mean the secondary dissociation of helium, the influence of the magnetic fleld on neutrino 
losses, and the spectrum of suprathermal emission of a photosphere in a wide energy 
range. 

For simplicity, we assume that the flreball is isotropic. However, all the results pre- 
sented remain valid also for jets with opening angles and Lorentz factors satisfying the 
condition 9 > 1/F. In the case of a jet, the source luminosity, mass-loss rate, etc., mean 
the corresponding values per unit solid angle multiplied by in. 



2 Neutrons in accretion disks 

Fission of helium nuclei due to inelastic collisions with high-energy ions in an accretion 
disk or with neutrons penetrating it is the main source of neutrons in accretion disks. It 
is usually assumed that the created neutrons are advected by the bulk plasma motion 
toward the central compact object^ with a velocity equal to the radial disk velocity [4, 
5] . However, if the accretion-disk plasma is tenuous enough, dynamical decoupling of the 
neutron and ion components of the disk takes place, which results in accumulation of 
neutrons. The term "decoupling" means the appearance of a non-negligible difference in 
the hydrodynamical velocities of the two components. 



2.1 Pile- up of free neutrons in a disk 

Consider the radial hydrodynamical motion of neutrons caused by the angular momen- 
tum loss due to neutron-ion (mostly, neutron-proton) collisions inside an accretion disk. 
Straightforward calculation yields 




where Vn_ and V^^^ are the orbital and radial velocities of the neutron component, Rg 
is the Schwarzschild radius of the black hole, and c is the speed of light. In the last 
expression, the time derivative of Vn was replaced by —Upn (Ki — Ki)/2, where i/pn is the 
rate of neutron-proton collisions and Vd is the orbital velocity of the disk. 

*That is, toward a neutron star or a black hole. Hereafter, for definiteness, we will mean a black hole. 
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The orbital velocity of the neutron component exceeds that of the disk, so that the 
excess centrifugal force should be balanced by the friction force due to the radial- velocity 
difference: 

Yl^ = ^ {Vy^ - vP) . (2) 

Here vj^^ is the radial velocity of the disk and it! is the current radius. Assuming 
Vd Vo, where Vq is the Keplerian velocity, we obtain: 

y^' = -\{'^)\y^'-yP)- (3) 

The quantity i^pn-R/Vo is equal to the number of collisions of a neutron during a single 
orbital rotation divided by 27r. 

The collision rate can be expressed in terms of the optical depth for neutrons r^, 
which is equal to the average number of collisions of a neutron penetrating the accretion 
disk perpendicular to its plane with a velocity about the Keplerian velocity. Assuming 
that the thickness of the ion and neutron "disks" are and /in, respectively, and that 
the cross-section of neutron-proton collisions is inversely proportional to their relative 
velocities, we find z/p^ = {(TV)pnnp (h^/K) = (Vq/K) {{aV)^^^ /Vq) Uphd = TdVo/{2hn). 
Here Up and are the number densities of protons and neutrons, respectively, and {aV)pn 
is the product of the proton-neutron cross section and their relative velocity averaged over 
the distribution function. Finally, the expression for the radial velocity of neutrons takes 
the following form : 

4+[i?/(2MPr| " ■ ^ ' 

The neutron pile-up and the formation of a neutron halo in the inner part of an accretion 
disk begins if Vj^^ becomes much less than Vd'\ i.e., ^ ^h^/R. Correspondingly, the 
accretion rate should be less than 



{aV) 



pn 



where nip is the proton mass. Here we used the relation V^ '^ = a (h^/ R^ Vo commonly 
accepted for a-disk models, while hd/R plays the role of the model parameter. If a = 0.1 
and hd ~ R, the limiting accretion rate Mmax corresponds to a luminosity about the 
Eddington one. 

As the accretion rate decreases, the role of the effect of neutron pile-up becomes more 
important, and eventually the number density of neutrons exceeds that of protons and 
the proper viscosity of the neutron component becomes important. Depending on the 
adopted disk modeL a further decrease in the accretion rate results in either stabilization 
of the ratio K^''VKi o^" the sign reversal of V^^^ and subsequent tending of its absolute 
value to zero. 
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2.2 Virialization of neutrons in the halo 



A neutron orbit is subject to precession caused by the gravimagnetic force by the rotating 
black hole. The precession is an effective mechanism of additional heating of neutrons 
and can result in the complete virialization of the neutron halo, which is possible if the 
angular velocity of the precession exceeds the rate of neutron-proton collisions, i.e., 

ac 



27ri?g \RJ^2 

where a is the angular momentum of the black hole in units of MRgc/2. Substituting the 
expression for i/pn in which the proton number density Up is given in terms of M and R, 
one can easily obtain the criterion for the accretion rate: 

M^^l^f^M^... (7) 



V27T V R 

It is clear that the considered mechanism can hardly provide virialization near the outer 
edge of the neutron halo even in the case of a rapidly rotating black hole (a — > 1). 
Nevertheless, virialization is possible in the inner part of the halo at distance 

/ 7.V \ 2/3 

2.3 Long-term variability 

Neutron-halo virialization near the black hole results in the formation of an almost 
isotropic wind that carries away about one-half of the neutrons created due to helium 
fission in the accretion disk and moves at a velocity comparable to the speed of light. 
Neutrons in the wind decay, and the resulting protons exert dynamical pressure on the 
accreted matter. At a distance less than tnVn, where is the wind velocity and is the 
lifetime of a neutron, this pressure is equal to 

P.-^^, (9) 

where r/He is the mass abundance of helium. In the case of quasi-spherical accretion (e.g., 
from a stellar wind from a massive companion star) , the neutron wind can terminate the 
accretion if Pyj exceeds the dynamical pressure of the infalling matter at the sonic point 
R^ ^ (c/Voo)^ Rg, where V^o is the velocity of the stellar wind at a large distance from the 
black hole. The corresponding condition can be presented in the following form: 
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As soon as the condition ([lOD is satisfied, the accretion terminates, but the source remains 
active due to the matter accumulated in the disk. This matter will be exhausted during a 
time interval of the order of the characteristic time scale of angular-momentum transfer at 
the outer edge of the disk. Then the source switches off, and the next cycle of the accretion 
from the stellar wind begins. The quiescent period lasts until the matter from the outer 
edge of the disk reaches the region of neutron-halo formation, i.e., it is approximately 
equal to the active period. The above- described sequence of events results in the long- 
term variability of a period very sensitive to the size of the outer edge of the disk, i.e., 
to the specific angular momentum of the accreted matter. In some cases, this period can 
amount to a few hundred or thousand years. 



2.4 Hard gamma-ray emission from the neutron halo 

If the neutron halo is virialized, pions can be created at distance R ~ 3Rg due to neutron- 
proton and neutron-neutron collisions. Assuming that the average rate of pion creation 
= (cTtt^) averaged over velocity distributions of the nucleons is the same for both types 
of collisions, we find the total "pion luminosity" 

K ~ ( Y + /^-'^P 1 ^'En, (11) 

where i?^ is the distance at which pion creation becomes possible and E.,^ ^ 140 MeV is 
the pion rest energy. If helium nuclei dissociate completely, the neutron number density 
can be determined using the balance condition 
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where the loss rate of neutron ^loss is determined mainly by their scattering into the 
black hole. Determining from Eq. (|12D and Hp from the continuity equation and 
assuming hj^ = R, we obtain the "pion luminosities" of neutron-neutron and proton- 
neutron collisions: 



^loss 2T7ip 



(13) 



The limiting values of these luminosities 

«7-lo38ar7HeMio(^l erg/s, (15) 



'/id' ^ 



R 



^max ^ ^qs^m^^ , _^ 1 erg/s (16) 
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can be obtained assuming M — Mmax- Here Mio is the mass of the black hole in units 
of IOMq, where Mq is the solar mass. About one-half of the pion energy is transformed 
to 7-quanta and electron-positron {e~e^) pairs. This energy fraction is reprocessed by 
virtue of comptonization of soft X-ray radiation and the electromagnetic cascade. 

Electrons and positrons cool rapidly due to interaction with the X-ray radiation of the 
disk, and the equilibrium distribution of e~e^ pairs is characterized by a mean energy of 
about 50—500 keV. The spectrum of photons with energies above a few tens of keV (i.e., 
above the cutoff of the blackbody emission from the disk) has a typical comptonization 
profile. The electromagnetic cascade, if it develops, leads to the formation of a hard power- 
law "tail" with an index of about 2 (we mean the photon number per unit frequency range) . 
Both spectral features are frequently observed during hard states of accreting sources — 
black hole candidates. Moreover, one should expect hard radiation with photon energies 
up to 70 MeV which, however, is hardly detectable since it carries much less energy 
compared to the hard X-ray radiation. 

The neutron halo can also be an efficient source of radiation in the deuterium line 
at energy a; 2.2 MeV. Deuterium is created due to evaporation of neutrons from the 
halo and their subsequent radiative capture by protons in the cold outer disk. The total 
luminosity in the line is equal to j3rfH.eMEd/[2mp), where j3 is the fraction of neutrons 
entering the outer disk. This fraction is strongly dependent on the velocity distribution 
of the neutrons and the halo size. 

3 Neutrons in relativist ic flows 

In this section, we discuss the characteristics of the neutron component in relativistic (and 
super-Eddington) plasma flows from compact objects. We call an object compact if the 
characteristic dynamical time, which is approximately equal to 3-Rg/c, is small compared 
to the neutron lifetime. This condition excludes from our analysis supermassive black 
holes in active galactic nuclei if their mass exceeds 3 ■ IO'^Mq. Moreover, for the sources 
with moderate densities and temperatures at the wind base, e.g., microquasars, some seed 
neutrons are also necessary. A neutron halo can be a source of such seed neutrons. 

3.1 Conditions for the formation and existence of the neutron 
component 

Gamma-ray bursts are probably the only example of sources in which the physical con- 
ditions necessary for independent (spontaneous) appearance of neutrons in a relativistic 
wind are met. For a nucleon density in the wind base 10"^— 10^ g/cm^ typical for gamma- 
ray bursts, the dissociation temperature is about 0.7 MeV for helium and is even less for 
other nuclei. In fact, the temperature Tq in the source is about an order of magnitude 
higher and amounts to 3 — 10 MeV. According to the generally accepted point of view. 
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ultrarelativistic plasma flows giving rise to gamma-ray bursts are formed during either 
catastrophic events in neutron-star evolution or the collapse of central parts of supermas- 
sive stars [16-20] . Therefore, the matter in the base of the flreball is initially enriched by 
neutrons, and even ordinary dissociation of the nuclei composing this matter makes the 
ratio of number densities of neutrons and protons amount to n-a/n^ ^ 1. 

There is another mechanism leading to creation of free neutrons in gamma-ray burst 
sources, namely, mutual conversion of protons and neutrons due to the following weak- 
interaction processes: p -|- Pe ^ + e"*" and p + e" ^ n + 1/^- A time of about [21] 

is enough for thermal equilibrium to be established, in which the number densities of 
protons and neutrons are approximately equal. Here Gp ~ 10~^^ erg • cm^ is the Fermi 

constant, Q is the typical energy of interacting particles (Q ~ 3To for a thermal distri- 
bution function), and gp^ ^ 1.25. In the course of the wind acceleration, the temperature 
decreases inversely proportional to the radius, and the time scale for establishing equi- 
librium increases sharply. Quite soon, at a temperature exceeding (mn — rup )c^, the 
density ratio becomes independent of the temperature behavior and remains equal to 
rj = rin/up fa 1. 

When the plasma temperature falls to about 70 keV, the recombination of protons 
and neutrons becomes possible. The process comprises two stages. First, a proton 
and a neutron form a deuteron via the reaction p + n d + ^ and then the reac- 
tion chain splits in different channels leading to the synthesis of helium according to 
the resulting reaction 3d ^^He + p + n. The recombination rate is determined mainly 
by the rate of deuterium synthesis. The recombination is negligible under the condi- 
tion Tr = 2 {av)dnp/3 < c/{RrT^), where {av)^ ~ 5 • 10~^° cm^/s is the constant of 
the deuterium-synthesis rate, is the distance from the source at which recombination 
becomes possible, and Fj. is the Lorentz factor at this distance. 

The Lorentz factor of the fireball increases almost linearly, F ^ R/Rq, until its max- 
imum Ti = L/{Mc^) determined by the source luminosity L and the rate of baryonic 
mass loss M. The radius Rq in the fireball base is assumed to be about 3Rg. To calcu- 
late the recombination rate of protons and neutrons, we use the continuity equation that 
determines the dependence of the proton number density Up on the radius: 

_ L 

~ 47ri?2 (1 + T]) rremc^ ' ^^^^ 

where m is the nucleon mass. Since a temperature of 70 keV is reached in gamma-ray 
burst fireballs at the acceleration stage, it can be easily shown that the recombination 
becomes significant if the following condition is met: 

\ / Tr\ \ 1 , . 

> 1. (19) 



1 km; Vl.5 MeVy (1 + 77)F^ 
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If ?7 ~ 1, then the fraction of free neutrons left after the recombination is {t^ + 1)^^ as 
long asri. + l<r7/|l— r^l, while at larger values of this fraction becomes exponentially 
small (if 1] < 1) OT constant (if r] > 1). 



3.2 Decoupling of proton and neutron components 

When the neutron component in a fireball is considered, we focus on the case where 
neutrons at the base of the fireball collide with protons frequently enough to be advected 
by the bulk flow. In fact, this means that the time between collisions of neutrons with 
other nucleons is less than the acceleration time scale R/{cr), i.e., 

a^pVnpR>cr L ;^ 12 (1 + 77) T.^Ledd- (20) 

0"np 

The expression on the right-hand side is calculated at the base of the fireball where F ^ 1 
and R ~ 3Rg. Also, we take into account that the relative velocity of protons and 
neutrons and the cross-section of their collisions at the threshold of decoupling are given 
by f ~ c/2 and a^p ~ 6 • 10~^^ cm^, respectively. In the obtained formula, the luminosity 
is related to the Eddington limit Lgdd = 27rmc^_Rg/crT ~ 1-3 ■ 10^^{M/Mq) erg/s, where 
(Tt ~ 6.6 ■ 10~^^ cm^ is the Thomson cross-section. Note that hereafter the quantity L 
may not necessarily be equal to the actual luminosity of the object and may significantly 
exceed the latter quantity. 

It is clear that the neutrons are indeed advected by the plasma outflow and form a 
neutral component if the luminosity L is much greater than the Eddington luminosity. 
However, in the course of acceleration and expansion of the wind, the decrease in the 
neutron-proton collision rate[] is much more rapid than the increase in the acceleration 
timescale. As a result, upon reaching certain distance, neutrons move almost freely. If 
the fireball acceleration continues beyond this distance, the Lorentz factor of the proton 
component exceeds that of the neutron component, i.e., decoupling takes place. The 
criterion of the decoupling is identical to Eq. (pO]) calculated at the distance Rg ~ TiRq 



from the source, where the fireball acceleration is terminated within the framework of the 
one-fluid approximation. The obtained condition can be presented in a more illustrative 
form: 

;^ r, = n/n, n ~ 5 kev, (21) 

where Tq is the temperature at the base of the fireball. 

For a given luminosity, the distance at which the decoupling takes place decreases 
with decreasing M. Therefore, the terminal Lorentz factor Fn of the neutron component 
as a function of F^ reaches the maximum at F^ ^ T^, and then decreases as F^^''^. The 
maximum value is approximately equal to F*. The limiting Lorentz factor Fp of the proton 



^To unify terminology, we write "proton" though, bearing in mind that hchum is created due to the 
recombination, the term "ion" should be used in some cases hereafter. 
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component increases monotonically with increasing Ti and, if neutrino losses are negligible 
(see below), can be calculated using the energy conservation law: Fp = + 77 (F^ — Fn). 
The value of Fn can be approximately calculated using the formula 

r r4/3 
(Ft + 2.37F4) 

In the case of gamma-ray bursts, the value of the decoupling parameter Fp/Fn ranges 
typically from 1 to 10, which is much less than the limiting Lorentz factor of the proton 
flow. On the contrary, the decoupling parameter is close to Fp in the case of jets in 
microquasars, so that the neutrons remain nonrelativistic. 



(22) 



3.3 Secondary dissociation of helium 

Consider the case where both the condition of efficient neutron recombination and the 
condition of the transition of the fireball to the two-fiow regime are satisfied. Then, prior 
to the decoupling of the proton and neutron components, secondary dissociation of helium 
becomes possible. In the accelerated reference frame comoving with the plasma flow, the 
neutron is subject to the inertial force mc^dT/dR directed toward the central source. 
Under the action of this force, the neutrons acquire a certain regular velocity and the 
neutron component is heated due to collisions with protons and helium nuclei. 

Let i^np and i/^a be the rates of elastic collisions of neutrons with protons and alpha- 
particles, respectively. Coulomb collisions with electrons cool ions quickly to energies 
of the order of the photon temperature T, so that ions can be considered motionless 
compared to neutrons. Therefore, steady-state values of the regular velocity v and r.m.s. 
velocity v of neutrons satisfy the following equations: 



,dr 



— - H V 

2 25 



^ 2dF 
dR 



(23) 



Just before a collision, the radial velocity of a typical neutron is v + Sv/2, where Sv 



c2 (dF/di?)/(z/, 



np 



l^n 



is the velocity increment between collisions. Now it is easy to find 



the mean energy of collisions of neutrons with helium nuclei, which is given by 



2m 



V + v5v + 



{Sv) 



2 (z/np + J^n 



np 



+ 



''np 4z/nQ 

2~ ~25 



+ 



^np ~l~ ^nc 

au: 



'np 



+ 



m {5vy 



- J 

(24) 

This value is equal to 4.1 m {SvY and 6.85 m (Sv)'^ in the limiting cases of hydrogen and 
helium plasma, respectively. 
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If E^a exceeds the threshold of the reaction n + ^He— d + ^B., which is about 18 MeV, 
multiphcation of neutrons due to hehum fission becomes possible. [| In this case, the 
total number of neutrons increases by a factor of 3^, where N is the number of inelastic 
collisions per neutron which occur during a time of about the fireball-expansion time scale 
R/{cT). Since the total number of colhsions during this time is about Nt = c/{6v), we 
arrive at the estimate N = VdmsNt/ {v^^ + z/na) -C A^t- 

Consider the case where r] > 1 and the recombination proceeded almost completely, 
i.e., the fireball at the onset of the secondary dissociation comprises only neutrons and 
helium nuclei. In this case, A^t ~ 18, and the rate of collisions resulting in the dissociation 
is equal to paiss ~ O.li/na, so that the number of neutrons increases by a factor of 10. 
Of course, rigorous calculation should account for the fact that protons are created due 
to the dissociation, and the ratio N/N^, as well as the value of A^t, are not constants 
(both these quantities decrease). Due to the appearance of protons in the course of the 
secondary dissociation, about one-half of the neutrons remains bound in helium nuclei 
after the recombination. But if t] is significantly less than unity but still Tj. ^ 1, secondary 
dissociation is almost impossible. 

The effect of secondary dissociation is significant for fireballs in GRB sources. At the 
same time, the influence of this effect on the jets in microquasars is negligible since the 
plasma in these jets is predominantly hydrogen. We should specially mention the so-called 
soft gamma repeaters (SGRs), i.e., neutron stars which produce sporadic bursts of soft 
gamma-ray emission of duration about 0.1 s and energy 10^^ — 10^^ ergs. No plausible 
theory of this phenomenon has been proposed up to now. Quite often, however, SGRs 
are assumed to be scaled-down analogs of GRBs. In these objects, the temperature is too 
low to cause the dissociation of helium at the base of the fireball. At the same time, low 
hydrogen abundance should lead to secondary dissociation. Thus, the formation of the 
neutron component in SGRs is entirely determined by the rate of nonthermal processes 
producing the seed neutrons. One cannot exclude the possibility that free neutrons play 
a noticeable role in SGRs, too. 

3.4 Radiative processes related to decoupling 

The collisions of two nucleons may be inelastic if their kinetic energy exceeds the pion- 
creation threshold of about 140 MeV. After the decoupling, the proton (ion) component 
remains cold, while the neutron one is heated up to subrelativistic temperatures. More- 
over, the relative velocity of protons and neutrons becomes relativistic, as well. If neutrons 
dominate the composition of nucleonic matter at the onset of the decoupling, the number 
densities of neutrons and protons become roughly equal, since inelastic collision of two 
neutrons leads to the conversion of one of these particles to a proton with a probability 
close to unity. 

■'■In fact, helium dissociates completely, *He ^ 2n + 2p, since all other reactions have lo-wer energy 
thresholds and greater cross-sections. 
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As the collision energy increases, the total cross-section of inelastic processes increases 
rapidly and becomes comparable to the cross-section of elastic scattering, so that a nucleon 
on average undergoes one inelastic collision (one should exclude the case 77 1 because 
of the above-mentioned reasons). The probabilities of creation of neutral or charged pions 
in proton- neutron collisions are almost equal.| The decay of charged pions, tt'^^tt^) 
ix^{li~)+i'^{v^) — > e+(e~)+z/e(z/e)+t';x+t'^, results in the appearance of three neutrinos and 
a positron (electron), and each particle carries approximately the same energy of about 
35 MeV. Neutral pions decay into a pair of gamma-quanta of energies about 70 MeV. 

Therefore, the decoupling leads to neutrino emission, originating mainly in a region 
near the saturation radius i?s, and to the electromagnetic cascade which generates hard 
gamma-rays escaping from the photosphere. The distance from the source to the pho- 
tosphere, which is located at different radii for photons of different energies, is generally 
much greater than the saturation radius, so that the electromagnetic emission is much 
weaker than the neutrino one. Nevertheless, it is easier to detect the electromagnetic sig- 
nal, while detection of neutrinos is virtually impossible even in the case of such powerful 
sources as gamma-ray bursts. 

Emission of high-energy particles due to pion decay results in the development of an 
electromagnetic cascade in the fireball plasma [|1^]. When analyzing the cascade, neutral 
and charged pions can be treated in a unified way. After creation of two e~e^ pairs by 
photons generated due to the decay of a neutral pion, these species differ only by the 
amount of energy transferred to electrons and positrons. Namely, the photons generated 
by a neutral-pion decay create two e~e^ pairs, and, since the energies of the daughter 
particles are roughly equal, they are almost indistinguishable from the electrons and 
positrons born in charged-pion decays. 

In the case of gamma-ray bursts, the problem is further simplified due to the fact 
that the cascade develops in the one-step regime near the photosphere. This means that 
thermal photons, upon scattering (comptonization) on initial ultrarelativistic particles, 
gain enough energy to create secondary pairs that cannot further transfer to photons the 
energy necessary for creation of new pairs, and thus the cascade terminates. 

What happens to a photon after the comptonization depends on its energy. If this 
energy exceeds the cascade threshold e^, the photon is absorbed by another photon of 
smaller energy, and an e~e^ pair is created. Otherwise, scattering by thermal electrons, 
in which the photon energy decreases gradually, is the dominant process. The quantity 

is universal in the approximation of a one-step cascade and is about 3 MeV. The 
newly created electrons and positrons are decelerated to almost thermal velocities prior 
to annihilation. As a result, all their energy is transferred to the background blackbody 
radiation, which plays the role of a heat bath with a heat capacity exceeding that of the 

§ Among other inelastic reactions, processes leading to deuterium formation (e.g., p + n — > d + Tr") take 
place. The created deuterons have a good chance to avoid new collisions and, therefore, dissociation. As 
a result, the fireball becomes an efficient source of deuterium: up to 20% of the ejected nucleon matter 
is converted to deuterium. 
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fireball plasma. 



3.5 Neutrino losses 

Neutrino emission is interesting mostly due to the fact that it can consume a large fraction 
of the fireball energy. The neutrino losses increase if the pion momenta are changed 
partially or completely before the decay. This is possible, e.g., due to interaction with 
the magnetic field. The decay time of a pion, equal to 2.6 ■ 10~^ s, exceeds the inverse 
gyrofrequency if the magnetic field is greater than 600 G. The energy density of such a 
field is many orders of magnitude smaller than the thermal energy density of the plasma 
in the decoupling region. Since the requirements to the uniformity of the magnetic field 
are not stringent in this case, we suppose that this effect is always present. Let us describe 
the effect of the magnetic field by the parameter ^, the angle between the field lines and 
radius. This parameter varies from to 7r/2 for different models of gamma-ray burst 
sources and jets. 

Let us estimate the neutrino losses in the case of large decoupling parameter. If 
r » Fn, all the three types of pions are created with equal probabihties, and their total 
energy in the ccntcr-of-mass frame of the colliding nucleons is approximately equal to 



mc^ yr/Fji. Taking into account the effect of the magnetic field, we find that neutrinos 
in the laboratory frame carry away the energy 

„ , (i- Ji-4rn/rcos2^) 
fti mc^r^r-^^ ^ ^ ^ . 

The total energy of the emitted neutrinos per proton is equal to 

dR 



J e^o-crin— , (25) 



where o"c ~ 2 ■ 10^^^ cm^ is the cross-section of inelastic proton-neutron collisions. We 
divide the integration region into three parts: (i) from the decoupling radius to the 



point where sin^ = ^^2T^/T; (ii) from the previous point to the saturation radius Rg — 
FpRd/r^; and (iii) from Rg to infinity. With accuracy sufficient for an estimate, we adopt 
rin — Ty^pF^/ (2r^), where F = [R/R^ F^ in regions (i) and (ii), and F = Fp in region (ni). 
If sin^ < Y/2Fn/Fp, then the formula ~ mc^F/2 is valid in any region. In the opposite 

case, ~ (mc^ F^ F^^" -"^ sin^ ^^)/4 in regions (ii) and (iii). 

The result of integration divided by the terminal proton energy V^mc^ shows how 
large are the neutrino losses compared to the energy required to accelerate the fireball: 

^" ^In ^ + 0.7) , sin e < ^2Fn/Fp; 

f / 2 \ (26) 

sin2^ + -^ ln-^-0.7 , sin^ > ^2F,/Fp. 
i p V sm ' V f 



SE, ^ 0.15 (^^^ r)r^ x < 



F V F 
1 p \ i 
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Here Tn ~ 4 is the optical depth of the decouphng region for neutrons and the factor 
(rd/Fn)^ is equal to about 2.5 rather than to a value of 3'^/^ corresponding to the decou- 
pling point. This is explained by the fact that the neutron flow continues to accelerate 
after the decoupling. Because of the adopted approximations, Eq. ( p6D is only applicable 
if Tp/Fn ^ 4. Note that if r] ^ 1 and the magnetic field is perpendicular to the radius, 
then over 50% of the source power is lost due to neutrino emission. This fact should be 
taken into account in calculations. The energy of a neutrino for a typical source with 
Fp ~ 10^ and Fp/Fn ~ 4 is about 30 GeV in the limiting case of a small angle 6 and about 
200 sin^ 9 GeV in the opposite case. 

3.6 Spectrum of the electromagnetic cascade 

Comptonization of the thermal radiation near the photosphere proceeds in the classical 
regime. Thus, if scattering of photons in the plasma is not taken into account, the flux 
of the comptonized photons per unit energy has the power-law form 

F, ^ ^.2N^m^c^T~^/^e-^''^, (27) 

and cuts off at the energy ^max ~ 47fT. Here A^e is the injection rate of primary 
electrons and positrons,[] 7i is their Lorentz factor, and me is the electron mass. The 
scattering in thermal plasma results in the cutoff of the initial spectrum at energies 
above meC^/r, where r = cxT^e-R/F is the Thomson optical depth at a given radius and 
He ~ 0.4 (Fp/Fn) ["rrieC^ / {2T)Y/^ is the total number density of electrons and positrons 
determined by the balance between the production and annihilation of e~e+ pairs. Strictly 
speaking, if r < et/(2meC^), then higher-energy photons in the narrow transparency in- 
terval with the lower boundary near Ei/t can survive, but these photons are insignificant 
for the determination of the cascade spectrum at energies below mgC^. 

The radiation emitted at the maximum optical depth for a given photon energy dom- 
inates the low-energy part of the electromagnetic-cascade spectrum (e -C mgC^). Such 
photons are subject to additional energy losses due to adiabatic expansion of the plasma. 
In this case, the upper boundary of the spectrum varies in the same way as the temper- 
ature of the background radiation and decreases by a factor of (Rt/RY^^, where i?T is 
the radius of the Thomson photosphere. Thus, the spectrum at energies e -C mgC^ is 
dominated by the radiation emitted at i? ~ [e/{TneC^)Y^^ R^. This radiation corresponds 
to a photon flux per unit energy that is (Rt/RY^^ = {rrieC^ / e^^^ times greater than the 
contribution of photons emitted from the Thomson photosphere. This yields the spectrum 

oc e"^ with the index q = —31/18. 

The spectrum above meC^ observed at infinity is determined by the local radiation 
of the electromagnetic cascade at the (energy-dependent) photosphere level. We point 
out three energy ranges. The opacity at energies from rrieC^ to et is determined by the 

^The contribution of secondary e~e^ pairs becomes important at energies less than {et/rricC^) T. 
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Klein-Nishina scattering by electrons and positrons, and the photospheric radius is equal 
to -Rph ~ {rricC^ / e)^/^ Rt;. Since -Rph Eq. (p7|), we arrive at the spectral index 

q = —7/10. If the photon energy exceeds e ~ et, the two-photon absorption becomes dom- 
inant, and the photosphere shifts away from the source as Rph ~ {meC^ /et){£ /etY^^^^ Rt 
with increasing e. Here the exponent is determined by the spectral index in the region 
6 < rjieC?. lil < e/et < {et/rrieC^Y^/^^ ~ 10, the distance from the source to the photo- 
sphere is less than /2t; so that the radial variation in the temperature should be taken into 
account, whereas i?ph > -Rt and T = const at higher photon energies. Correspondingly, 
the spectral index q = —107/54 ?s —2 in the first case and q = —20/9 in the second one. 
The results obtained may be summarized as follows: 



' Q = 


-31/18, 


e < rrieC^; 


Q = 


-7/10, 


rricC^ < e < Et, 


Q = 


-107/54, 


et<e < min{10£:t,£:max}; 


. Q = 


-20/9, 


min{10£t,emax} <e < ^max 



(2^ 



The main energy- carrying part of the spectrum (pSD is the region between St and £up = 
min{10 et, ^max}- Multiplying by Fp, we obtain an energy range of 2 to 50 GeV in the 
laboratory frame. 

Applying the same method that was used to obtain Eq. (^), we calculate the fraction 
of energy 5-Ecas corresponding to the emission of the electromagnetic cascade. In this 
case, the integration is performed from i?T to infinity, the effective value of sin^ 6 is equal 
to 1, and the obtained result is multiplied by the ratio of the total energy of photons 
in the considered range to the energy of electrons and positrons injected at the distance 



i?T- This ratio is about \n{£u^/£t) y^up/^max- For large decoupling parameter we have 
f^max > 10 et, which yields 

The effective area of the gamma-ray telescope onboard the GLAST satellite will allow 
one to observe the electromagnetic-cascade emission from gamma-ray bursts at redshifts 
z 0.1 provided the burst energy is 10^^ erg. 

The energy is partially emitted in the form of gamma-ray quanta directly coming 
from the decay of neutral pions. These photons fall in the 100-GeV energy range in the 
laboratory frame. The fraction of energy carried by such a radiation is given by the 



following formula obtained in |15]: 



mcrxFn^e'^e Ve 

Here Si is the average energy (in the fireball frame) of gamma-quanta born in the de- 
cay of neutral pions. This formula is corrected for the more detailed model of the 
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electromagnetic-cascade spectrum developed above. However, we find that the numer- 
ical value of 6Ej is almost unaffected by the specific form of the cascade spectrum. This 
value is changed by a factor less than 1.3 compared to the simplified model. The photons 
from ttq decay can be detected by ground-based Cerenkov telescopes in the near future 
when such facilities as MAGIC or ARGO become fully operational. 

4 Conclusions 

It is shown above that, in certain cases, free neutrons are inevitably created in astrophysi- 
cal sources of X- and gamma-ray radiation due mainly to dissociation of helium nuclei. A 
temperature of 0.7 MeV required for this process is either achieved in the source initially 
(gamma-ray bursts, inner regions of accretion disks in advection-dominated models) or 
reached due to the large relative velocity of seed neutrons and ions (e.g., jets in micro- 
quasars and SGRs). 

A neutron halo in the inner region of an accretion disk originates owing to the effect 
of neutron pile-up. In this case, the radial velocity of the neutron component is usually 
much less compared to the electron-ion plasma, while the temperature of this component is 
high enough to ensure coUisional dissociation of helium nuclei. Neutron- neutron collisions 
give rise to the neutron wind whose dynamical pressure is enough to terminate accretion 
from the wind of the companion star. As the accretion is terminated, the neutron wind 
weakens, the accretion starts again, and the entire automodulation cycle of the accretion 
rate is repeated. A fraction of neutrons from the wind is intercepted by the cold outer 
disk where these neutrons lead to the formation of nuclear lines, primarily the 2.2-MeV 
deuterium line. In the immediate vicinity of the central object, the neutrons undergo 
inelastic collisions and produce pions, thus initiating the electromagnetic cascade which 
finally produces the comptonized radiation with a hard power-law spectrum and the 
annihilation line at 511 keV generated by positron annihilation. 

In the course of acceleration of a plasma flow comprising a neutron component, the 
neutrons accelerate due only to their collisions with protons. As soon as the time between 
collisions exceeds the characteristic time scale of the fireball expansion, the acceleration 
of neutrons terminates. Thus, the Lorentz factors of the neutron and proton components 
can become strongly different. Such a decoupling results in the dissociation of helium 
nuclei (if helium exists in the flow) due to interaction with seed neutrons and leads to 
an exponential increase in the number of neutrons. Moreover, the decoupling gives rise 
to inelastic nucleon collisions in which about one pion per neutron is born. About a 
half of the pion energy is lost for neutrino emission, which has an important effect on 
the energy balance of the source. Pion decay gives rise to an electromagnetic cascade in 
which a signiflcant amount of electron-positron plasma is created. Due to the cascade, the 
fireball photosphere shifts outward, and hard gamma-ray emission appears. Observations 
of the spectral features of the electromagnetic cascade allow one to determine the plasma 
parameters and can be indirect evidence for the existence of neutrons in the fireball. 
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The above-mentioned effects are also possible in the absence of decoupling if the fire- 
ball is strongly inhomogeneous (variability time of the source about i?g/c) or strongly 
anisotropic (with a scale about l/F) and its Lorentz factor exceeds F*^^. In this case, 
the neutrino losses decrease, while the electromagnetic-cascade spectrum remains almost 

unchanged. 

Finally, since neutrons carry a significant fraction of the fireball energy which is con- 
verted to radiation only upon neutron decay, they can strongly affect the light curves 
of transient sources, including the formation of a secondary maximum and a bimodal 
distribution over the duration [14, 15]. 
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